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We propose and experimentally verify a scheme to engineer arbitrary states of traveling light 
field up to the two-photon level. The desired state is remotely prepared in the signal channel of 
spontaneous parametric down-conversion by means of conditional measurements on the idler channel. 
The measurement consists of bringing the idler field into interference with two ancilla coherent states, 
followed by two single-photon detectors, which, in coincidence, herald the preparation event. By 
varying the amplitudes and phases of the ancillae, we can prepare any arbitrary superposition of 
zero- one- and two-photon states. 



PACS numbers: 

Modern quantum information science relies upon two 
key concepts - superposition and entanglement. As fun- 
damental tenets of quantum mechanics, they govern the 
way information can be shared, transferred or measured. 
This information lives in a Hilbert space, in the form 
of a quantum state. Consequently, the ability to ac- 
tively control the coherent dynamics of a quantum state 
is paramount for quantum information technology. This 
task forms the essence of quantum state engineering 
(QSE). In the optical domain, a widely-used approach 
for QSE involves generating a "primitive" quantum state 
and then manipulating it, for example, by bringing it 
into interaction with an ancillary system. Employing ap- 
propriate measurements, the ancilla is then traced out, 
leading to reduction of the overall system to the desired 
target state, ready to be detected and characterized. 

In modern quantum optics, the "primitive" is com- 
monly the state of correlated photon pairs produced in 
spontaneous parametric down conversion (SPDC). Con- 
ditional photon detection on one or both channels is then 
employed to produce the state of interest. This technique 
has been successfully applied to engineer complex entan- 
gled states of dual-rail optical qubits [1] , albeit mostly in 
a postselected manner: we do not know that the state 
has been prepared until it is detected and destroyed. A 
variety of single-mode states have been prepared without 
resorting to postselection: photon number (Fock) states 
|1) [2] and |2) [3], single-rail optical qubits [4, 5 , photon- 
added states [6 and the "Schrodinger kitten" state [7]. 
However, engineering of arbitrary quantum states of light 
has not yet been demonstrated. 

A fundamental impediment to optical QSE arises be- 
cause equidistant energy levels of a harmonic oscillator 
(such as an electromagnetic mode) cannot be individu- 
ally accessed using classical control signals. Outside the 
optical domain, this challenge has been addressed in ion 
traps [8] and high-Q microwave cavities [9 by coupling 
the oscillator energy eigenstates to a two-level atomic or 
spin system [ID] . Very recently, arbitrary superposition 
of Fock states were synthesized inside a superconduct- 



ing cavity by means of coupling to a Josephson phase 
qubit [11 . However, the loss of coherence in this scheme 
provides a fundamental limitation to the obtainable ac- 
curacy and complexity of the prepared state. 

While traveling field implementations do not suffer 
from this kind of decoherence and also automatically sat- 
isfy DiVincenzo's criterion of "flying qubits" [12], they 
are still beset with the problem of inefficiencies and losses. 
There exist a number of theoretical proposals for imple- 
menting optical QSE [reviewed in detail in [13]], for ex- 
ample, using coherent displacements and photon subtrac- 
tion operations [14] |15] [16] , repeated parametric down- 
conversion [17] , or continuous- variable postselection [18] . 
Here we report, for the first time, a postselection- free 
experiment on synthesis of single-mode coherent super- 
positions of Fock states up to the two-photon level with 
stable phase relations and a high overall efficiency. 

remote state preparation 




FIG. 1: Scheme to produce arbitrary quantum-optical states 
at the two photon level. The idler output of parametric down- 
conversion is mixed with two weak coherent states on beam- 
splitters, the coincidence detection from the two detectors 
heralds the production of the expected state in the signal 
channel. SPCM, single-photon counting module. 

Figure 1 shows a simplified layout of our experiment. 
The two-mode state produced by down-conversion can 
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be described (neglecting normalization) as 

|*) = |0> + 7 |1„ h) + 7 2 \2„ 2i) + 0(7 3 ), (1) 

where 7 2 is the probability of generating a single photon 
pair and the signal and idler modes are represented by 
s and z, respectively. The conditional measurement is 
performed by mixing the idler photons with two weak 
ancillary coherent states on symmetric beam splitters, 
and detecting photon arrivals in two output channels. 
When the two detectors "click" simultaneously, the non- 
locally prepared signal state is a superposition of |0) , |1) 
and |2). 

A qualitative explanation is as follows: each of the two 
trigger photons can originate either from SPDC or from 
the ancillary states. Since SPDC emits the same num- 
ber of photons into both channels, a coincidence detec- 
tion event indicates that the photon number in the signal 
path is 0, 1 or 2. Furthermore, if the optical modes of 
the ancillary and idler fields are matched spatially and 
temporally, the origin of the photons giving rise to reg- 
istration events is fundamentally indeterminable. As a 
result, the prepared state p is not simply a statistical 
mixture, but a coherent superposition of the first three 
Fock states. The coefficients of this superposition can be 
controlled with the amplitude and relative phase of the 
two ancillae. 

The commercial detectors employed in state prepa- 
ration (Perkin-Elmer SPCM-AQR-14-FC) "click" in re- 
sponse to single photons, but are unable to resolve their 
exact number. However, the parameters a, /?, 7 are on 
the scale of 0.1 [19], so the probability that a given "click" 
has occurred in response to two or more photons enter- 
ing the detector is two orders of magnitude lower than 
that due to a single photon. Thus it is safe to assume 
that all "clicks" are associated with single photons, and 
the contribution of n > 3 terms in the signal channel is 
insignificant. 

It is interesting that the entire state preparation occurs 
remotely, without any manipulation of the signal channel 
[20] . This is possible because the SPDC output is 
entangled. The low degree of entanglement associated 
with small 7 does not restrict the range of quantum states 
that can be engineered, but only reduces the frequency 
of successful preparation events. 

In the limit of small a, (3 and 7, the state prepared in 
the signal channel can be calculated to be 

|^) = ao|0)+ai|l) + 02|2) (2) 

with 

a 2 a/3 /?7 . . j 2 

A close examination of Eq. ([2| reveals an interesting fea- 
ture: complete elimination of the |0) or |1) components 
is possible simply by blocking one of the ancilla fields. If 



a = 0, the first SPCM will trigger only when the first 
photon comes from SPDC and thus the signal channel 
would always contain at least one photon. If we instead 
block the second coherent state, a coincidence event will 
occur if the two triggering photons come from the first 
beam splitter. Due to the Hong-Ou-Mandel effect [21], 
this cannot happen if one photon originates from SPDC 
and one from the coherent state \a). Hence a coinci- 
dence event with (3 = cannot occur if exactly one pair 
was generated in SPDC. 

In a complete theoretical description, several param- 
eters must be taken into account in order to reproduce 
and validate our experimental results. First, remote state 
preparation using down- conversion leads to generic inef- 
ficiencies (due to optical losses, mismatch in the modes of 
the signal and the local oscillator used in detection, dark 
counts of detectors, inefficient photodiodes and electronic 
noise) [2], which can be treated as partial absorption of 
the signal light. Second, imperfect mode-matching of the 
ancilla fields with the idler photons leads to partial dis- 
tinguishability and loss of coherence among Fock terms 
in the prepared state. Finally, phase drifts and higher 
photon number components must be taken into account. 

We worked with 1.7-ps laser pulses at 790 nm emitted 
by a mode- locked Ti: Sapphire Coherent Mira 900 laser 
with a repetition rate of ~76 MHz. Most of the laser 
energy was directed towards a LBO crystal for frequency 
doubling. This provided us with a strong pump beam 
at 395 nm, which drove SPDC in a periodically poled 
KTiOP04 crystal. This crystal was phase-matched for 
type II down- conversion in which the two output modes 
were spatially and spectrally degenerate but had orthogo- 
nal polarizations. They were then separated spatially on 
a polarizing beam splitter. Characterization of the signal 
field was performed using homodyne tomography [22] , 
with the local oscillator (LO) obtained from the main 
output of the laser. 

The ancilla fields also originated from the master laser. 
As evident from Eq. stability of the relative phase 
arg(a) — arg(/3) of these states was critical for preserv- 
ing coherence among different number states. The re- 
quired stability was achieved by employing calcite beam 
displacers in an interferometerically stable configuration 
[23] . The relative phase was found to stray by less than 
0.04 radians on a typical time scale of a data acquisition 
run (around 30 minutes). 

On the other hand, the stability of the common phase 
[arg(a) + arg(/?)]/2 of the two ancilla states as well as the 
phase of the SPDC pump field arg(7) was not critically 
important. The only effect of the drift of these phases 
was to modify the optical phase of the signal state 
The latter was continuously monitored in the process of 
homodyne detection, as we describe below. 

The homodyne detector employed two Hamamatsu 
S5972 photodiodes (85% quantum efficiency) and fea- 
tured a 90-MHz bandwidth, which permitted time- 
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resolved quadrature measurements at the LO pulse rep- 
etition rate [19]. Time- integrated photocurrent samples 
associated with each LO pulse were normalized with re- 
spect to the vacuum state, yielding quadrature data of 
the state measured. For each state reconstruction, 50000 
samples were typically acquired. 

In order to reconstruct the state, we also needed to 
know the phase of the signal in reference to LO for each 
quadrature sample. This information could be deduced 
by analyzing time-dependent quadrature statistics, be- 
cause at any given phase <\> the average quadrature value 
behaves as (Q^) ex sin0. However, in our experiment 
the phase varied on a time scale of ~ 0.5 rad/s, so the 
coincidence-triggered data, acquired at a rate of 20-200 
Hz, were insufficient for reliable phase reconstruction. 
Therefore we additionally acquired the homodyne signal 
triggered by events from both individual single-photon 
detectors, which occurred at a much higher rate of 25- 
100 kHz. The signal states prepared by these events are 
phase-sensitive coherent superpositions of Fock states |0) 
and |1). Average quadrature values of these states over 
~ 0.06 s periods were used to determine the LO phase at 
each moment in time. The phase-quadrature pairs from 
the acquisition were then fed to a maximum-likelihood 
reconstruction algorithm [24] to estimate the density ma- 
trix of the signal state. 
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FIG. 2: Results for Fock states |0) (a), |1) (b), and \2) (c). 
The Wigner function and the density matrix (absolute values) 
are displayed. All state reconstructions feature correction for 
55% detection efficiency. 



Each acquisition run began with preparation and re- 
construction of single-photon Fock states, by blocking 
the ancilla fields (i.e. a = f3 = 0) and triggering the 
quadrature acquisition by individual SPCMs. This al- 
lowed us to test the performance of our setup and deter- 
mine its overall efficiency. The latter was typically found 
to be between 0.52 and 0.55, in agreement with the ex- 
pected value of 0.54 calculated as a cumulative quantum 
efficiency of all parts of the experiment [19]. In subse- 
quent applications of the state reconstruction algorithm, 
we corrected for this known inefficiency using a value of 
rj = 0.55 [24 . In order to further test our procedure, 
we have prepared and reconstructed the three basis Fock 
states (Fig. 2). 

Figures 3, 4, and 5 display coherent superpositions 
of states |0) and |1), |0) and |2), |1) and |2), respec- 
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FIG. 3: Reconstructed superpositions of states |0) and |1). 
The single-photon fraction increases from left to right. 
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FIG. 4: Reconstructed superpositions of states |0) and |2). 
The two-photon fraction increases from left to right. 



tively. The states in Fig. 3 were prepared by condi- 
tioning the quadrature acquisition on click events from 
only the first SPCM in the presence of a single ancilla 
field | a). With this setting, the heralded signal state 
was = a |0) + 7|1) [4]. The states in Figs. 4 and 
5, on the other hand, were prepared by blocking one 
of the weak ancilla inputs and triggering the homodyne 
acquisition on coincidence counts. The three examples 
in each figure contain varying weights of superposition 
terms, implementing gradual transition from one Fock 
state to another. It is interesting to note that the state 
in Fig. 4(a), i.e. the vacuum state with a small contribu- 
tion of 1 2), is a good (95% fidelity) approximation of the 
even "Schrodinger kitten" state |a) + |— a) with a = 0.60. 

Superpositions of all three Fock states are displayed 
in Fig. 6. In this acquisition run, the phases of ao, 
and CL2 were set equal by adjusting the classical interfe- 
rence of the two ancilla fields in the non-blocked output 
of the second beam splitter to a constructive fringe. In 
the three states displayed, the values of (3 and 7 were 
maintained approximately the same but the amplitude a 
of the first coherent state was varied. As a result, we ob- 
serve a gradual transition between the vacuum state and 
an equal- weight coherent superposition of |1) and |2). 

In all of the above examples, the zero phase reference 
of the reconstructed state could be chosen to make all 
density matrix elements real within the experimental er- 
ror, leading to reflection symmetry of the reconstructed 
Wigner functions. This is no longer the case in our two 
final acquisitions, aimed to demonstrate our control over 
the phases of (Fig. 7). Here the ancillae's phases were 
adjusted to observe interference midway between a bright 
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FIG. 5: Reconstructed superpositions of states |1) and |2) 
The two-photon fraction increases from left to right. 
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FIG. 7: Superpositions of states |0), |1) and |2). The mag- 
nitudes of ao,ai, and a2 in (a) are approximately equal to 
those in (b). Amplitudes ao and a\ are real whereas the a2 
are complex conjugate of each other in (a) and (b), so the 
Wigner functions appear as mirror reflections of each other. 
The real and imaginary components of the density matrix are 
shown for both states. 
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FIG. 6: Reconstructed equal-phase superpositions of states 
|0), |1) and 1 2). In all plots, a\ ~ a<2,\ the vacuum amplitude 
ao decreases from left to right. 



and dark fringe, on opposite fringe sides in Figs. 7(a) and 
(b). This implies a phase difference of ±7r/2 between a 
and (3. If the optical phase of the reconstructed state is 
chosen so that a\ and a2 are real, the amplitude of ao 
in Fig. 7(a) is complex conjugate of that in Fig. 7(b). 
As a result, the two Wigner functions themselves are not 
symmetric, but exhibit mirror symmetry with respect to 
each other. 

In order to estimate the quality of our state prepara- 
tion, we fitted each of the reconstructed states p with 
Eq. ([2|, aiming to maximize the fidelity F = (i/;\p\i/;). 
The fit parameters were restricted to agree with the ex- 
perimental conditions under which each data set was 
taken. The obtained best fit parameters are shown in 
Figs. 3-7 next to each plot. All states studied in this 
work exhibited at least a 76% fidelity with the fits. 

To summarize, we have experimentally demonstrated 
production and characterization of superpositions of Fock 
states |0), |1) and |2) using parametric down-conversion 
and non-local state preparation. We used two weak co- 
herent states as ancillae, which we mixed with the SPDC 
idler channel in a phase-stable setting to devise a mea- 
surement that remotely prepares the desired signal state. 
Blocking one of the coherent states allowed us to remove 
the |0) or |1) component from the superposition. Exten- 
sion of this method to a higher photon number domain is 
possible, however requires additional theoretical work as 
well as enhancement of experimentally accessible photon- 
pair sources. 
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